Raman imaging of doping domains in graphene on Si02 
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We present spatially resolved Raman images of the G and 2D lines of single-layer graphene flakes. 
The spatial fluctuations of G and 2D lines are correlated and are thus shown to be affiliated with 
local doping domains. We investigate the position of the 2D line - the most significant Raman peak 
to identify single-layer graphene - as a function of charging up to |n| ~ 4 x 10^^ cm~^. Contrary to 
the G line which exhibits a strong and symmetric stiffening with respect to electron and hole-doping, 
the 2D line shows a weak and slightly asymmetric stiffening for low doping. Additionally, the line 
width of the 2D line is, in contrast to the G line, doping-independent making this quantity a reliable 
measure for identifying single-layer graphene. 

PACS numbers: 63.20.K-, 78.30.-j, 78.30Na, 81.05.-t, 81.05.Uw 



Graphene has attracted increasing attention over the 
last few years [l|, [H . Its unique electronic properties 0, 0] , 
mainly due to the linear energy vs. momentum disper- 
sion and the electron hole symmetry near the charge 
neutrality point, makes it an interesting nanomaterial 
for high mobility electronics [H, 0] . Raman spectroscopy 
has proven to be a powerful tool to distinguish single- 
layer graphene from few-layer graphene and graphite 
B S HO, [HI- The particular electronic structure of 
graphite and graphene leads to Kohn- anomalies in the 
phonon dispersion at the F and K points [12, 13]. At F, 
the Born-Oppenheimer approximation, which is usually 
employed for the calculation of phonon frequencies, is 
no longer valid because of the Kohn-anomaly [14]. This 
leads to a pronounced stiffening of the Raman G line 
upon positive or neg ative charging (p or n-doping) of the 
graphene sheet [Hf [H, [13 • In this paper, we investi- 
gate how spatially resolved Raman spectroscopy can be 
used to probe doping domains and local charge fluctu- 
ations. While electron-hole puddles (i.e. local charge 
fluctuations) have been predicted to be responsible for 
the finite conductance at vanishing (average) charge car- 
rier density [l^ and have recently been observed using a 
scanning single electron transistor [19], the identification 
of different doping domains might be desirable to investi- 
gate novel graphene devices. Here, we report on Raman 
measurements on back gate induced charged graphene 
and on Raman imaging of doping fluctuations of isolated 
graphene flakes. We focus on the correlation between 
the shifts of the G line and the 2D line (or D* line). The 
latter one is the most significant Raman peak in single- 
layer graphene 0, 0. Within the low charging regime 
(up to ±4 X 10^^ cm~^) obtained in our experiments, the 
2D line stiffens for both electron and hole charging while 
its line width (in contrast to the G line) is not affected 



by charging. A good correlation between the shift of G 
and 2D lines is observed. However, the spectral resolu- 
tion and lateral resolution are not sufficient to resolve 
electron- hole puddles as shown in Ref. 0. Therefore 
we refer to charging (i.e. doping) domains rather than to 
local electron-hole puddles. 

We present both Raman images of isolated graphene 
and Raman measurements on electrically contacted 
single-layer graphene on 300 nm Si02, where highly 
doped Si is used as back gate. The samples are prepared 
by micromechanical cleavage [2Q|] and Raman imag- 
ing [21] is used to select single-layer graphene flakes 0,[3- 
By electron-beam lithography we pattern scanning force 
microscopy pre- mapped electrodes (5 nm Cr/60 nm Au) 
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FIG. 1: (color online) (a) Schematic illustration of a graphene 
sample with incident laser, (b) Raman image of a contacted 
single-layer graphene flake, where the integrated 2D line 
(green; dark) and the inverse Si line (at approx. 520 cm""*^), 
which is screened by metal electrodes (yellow; bright) are col- 
ored, (c) Raman spectra of the G and 2D line of charge- 
neutral and charged graphene measured at spot P marked in 
(b) at T = 295 K. 



on the graphene flakes, which flnahy ahows to apply a 
back gate voltage Vg between the Si++ substrate and 
graphene. The experimental setup is shown in Fig. 1(a). 
The Raman data are recorded by using a laser excita- 
tion of 532 nm (£^^=2.33 eV) through a single- mode op- 
tical flber, whose spot size is limited by diffraction. A 
long working distance focusing lens with numerical aper- 
ture of approx. 0.80 is used to obtain a spot size of 
approx. 400 nm. We use a laser power below 2 mW 
such that heating effects can be neglected [22|. A Ra- 
man image of a measured device is shown in Fig. 1(b) and 
the Raman spectra corresponding to point P (Fig. 1(b)) 
for a charged and charge-neutral case are plotted in 
Fig. 1(c). The G hue shifts approx. by 3 cm~^ due to hole 
charging (of n ^ —4 x 10^^ cm~^) and a corresponding 
(hardly visible) small shift in the 2D line is observed, too. 
Fig. 2 shows Raman data for varying back gate voltages 
Vg (Figs. 2(a,b)), which by utilizing a simple capacitor 
model can be substituted by the electron/hole concen- 
tration n = a{Vg-Vg^). Here a ^ 7.2 x 10^° cm-^/V g 
and 2.5 V marks the charge neutrality point, which 

has been determined by transport measurements. A typ- 
ical back gate characteristic of the investigated device is 
shown in Fig. 4 in Ref. [1^. The symmetric hole and elec- 
tron charging-dependent stiffening of the E2g F-phonon 
was recently explained as the effect of non-adiabaticity 
[3, [H, [H, [l3 5 the fact that the time-scale of phonon 
oscillations in graphite is not long compared to the elec- 
tron relaxation time (see e.g. Eq. (6) in Ref. 0], and 
dashed line in Fig. 2(c)). Results of time-dependent per- 
turbation theory were, however, so far only presented 
for phonons at the F point. For the phonon between K 
and M which is responsible for the 2D line, it has been 
argued Q that non-adiabatic effects are negligible and 
the influence of charging can be reproduced by a stan- 
dard adiabatic phonon calculation. While non-adiabatic 
calculations predict a phonon stiffening for hole-doping 
and a phonon softening for electron-doping [ll, [lH , we 
observe (see Fig. 2(d)) a slight (asymmetric) stiffening 
for both electron and hole-doping in agreement with the 
measurements of Ref. [17]. 

Another significant difference between the G and 2D 
lines is their line width (FWHM, shown in Fig. 2(e)). The 
G phonon {qc ^ 0) shows a rather strong change as func- 
tion of carrier density [JB^Il^] which is due to the fact that 
the Pauli exclusion principle prevents the phonon from 
decaying into an electron- hole pair for \Ef\ > Tiuog/'^i as 
illustrated in Fig. 2(f) [13] • The decay of the dispersive 
D phonon with large wave vector qd '^^ unaffected by the 
Pauli exclusion principle for low doping (Fig. 2(g)). It is 
expected that the 2D line width stays constant up to a 
charging that corresponds to a Fermi level shift as large as 
the exciting laser energy \Ef\ ^ El. This has a practical 
implication: Since the peak width of the 2D line, which 
has been recognized as the most striking feature to dis- 
tinguish single-layer from few-layer graphene 0, B Q , is 
insensitive to doping, it is a reliable - doping-independent 
- measure for identifying single-layer graphene. 
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FIG. 2: (color online) Raman spectroscopy of the G and 2D 
line of charged graphene. (a) and (b) show two-dimensional 
maps of the G and 2D line, respectively as function of ap- 
plied back gate voltages. Both plots are peak height normal- 
ized, (c) shows the G line peak shift [extracted from (a)] as 
function of the induced carrier concentration. Dashed lines 
show theoretical expectations for finite temperature [l3] . (d) 
same for the 2D peak shift. The same relative frequency- 
scale as in (c) is used, indicated by the horizontal scale bars 
(1 cm~^ /phonon). (e) shows the normalized band width 
(FWHM) of G and 2D as function of doping, (f) Schematic 
illustration of the presence (left panel) and absence (right 
panel) of a F phonon decay into an electron- hole pair, (g) 
Presence and absence of D-phonon decay into an electron- 
hole pair. Ki and K2 denote the two inequivalent iT-points. 



Spatially resolved Raman spectroscopy (Raman imag- 
ing) sensing the doping-dependent G and 2D line shifts 
provides an interesting tool to investigate charge fluctu- 
ations and doping domains in graphene. Here we present 
an example of Raman images (80x45 pixels) of a graphene 
flake on Si02 (Fig. 3), where most attention has been 
paid to the G and 2D peak positions and widths, their 
fluctuations and cross-correlations. The average G peak 
position (Fig. 3(c)), uoq-, measured on the flake shown 
in Fig. 3(a), is 1585.4 cm~^ and the root mean square 
(RMS) of the peak fluctuations is 3.3 cm~^. This fluc- 
tuation also nicely explains the wide spread range of uoq 
reported in the literature 0, i, i, [10, [HI. The 2D line 
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FIG. 3: (color online) Raman images (RIs) of an isolated 
single-layer graphene flake on Si02. (a) Green filtered opti- 
cal image (the two white dots are metal markers), (b) RI 
of the FWHM of the 2D line, proving the single-layer char- 
acter of the flake. The white dashed line marks the metal 
marker. (c,d) shows RIs of the G and 2D peak positions, 
respectively. Most attention is paid to the correlated peak 
fluctuations, which are due to doping fluctuations. The white 
arrows mark correlated low doping areas, whereas the dot- 
ted box highlights a region with low charging fluctuations, 
(e) [(f)] Two-dimensional cross-correlation of the fluctuations 
of the G and 2D line positions [widths]. In both cases the 
gray scale is adjusted by setting the maximum value to black. 
Strong correlation results (e) therefore in a centered peak at 
(0,0), whereas low concentration leads to a distributed func- 
tion [e.g. (f)]. 



(Fig. 3(d)) is centered around 0020 = 2679.4 cm~^ and 
fluctuates with 0.9 cm~^. The ratio of the fluctuations 
of the G and 2D line agrees well with the ratio of the 
doping-dependent G and 2D stiffening (Figs. 2(c,d)) of ~ 
3.2. In Figs. 3(c,d) doped regions on the imaged flake 
can be observed [11]. In the upper part we see that 
towards the edges of the graphene sample charging is 
suppressed, whereas in the entire inner area signifi- 
cant charging is present. Focusing on a quite uniform 
area (dotted box in Figs. 3(c,d)) we find that the uoq- 
fluctuations are approx. 0.6 cm~^, which corresponds 



to An ^ 2.4 X 10^^ cm~^. Please note, that the lo- 
cal charge fluctuations due to electron-hole puddles pre- 
sented in Ref. [l^ are by one order of magnitude smaller 
in amplitude and have been measured at low tempera- 
ture and controlled environment. In addition, our lateral 
resolution is limited by the laser spot size. Figs. 3(c,d), 
however, nicely illustrate the correlation of G and 2D 
fluctuations (see e.g. white arrows). This can be also 
expressed by the two-dimensional cross-correlation func- 
tion C^'^^(Ar) = j uocif -\- Ar) uj2D{r)dxdy, which ex- 
hibits a clear narrow central peak (Fig. 3(e)). In contrast 
^FWHAii^^) ^^^^ show no measurable correlation of the 
width of G and 2D (Fig. 3(f)) as expected, since they are 
uncorrelated quantities. Indeed the FWHM of the 2D 
line does not exhibit any doping dependence as shown in 
Fig. 3(b). The line width of the 2D line stays constant at 
33±0.08 cm~^ testifying its good quality for identifying 
single-layer graphene. 

In summary, we have presented Raman shifts of the 
G and 2D line for back gate induced charged graphene. 
We have discussed the spatial variations in the G and 
2D peak positions which can be attributed to different 
doping domains. In the low doping regime no clear dis- 
tinction between electron and hole-doping can be made 
since for both G and 2D line, stiffening is observed. How- 
ever, absolute doping fluctuations can be estimated. This 
technique is promising to investigate, e.g., electric field 
distributions in (side) gated graphene devices [26l |. 
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Note added — During completion of the manuscript, 
we became aware of a quite similar work [l^ in which 
doping-dependent fluctuations and a correlation between 
the position of G and 2D line are observed during 
non-spatially resolved measurements on many different 
graphene flakes. 
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